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Abstract 

We apply Tsallis's ^-indexed entropy to formulate a non-extensive random matrix theory (RMT), 
which may be suitable for systems with mixed regular-chaotic dynamics. The joint distribution 
of the matrix elements is given by folding the corresponding quantity in the conventional random 
matrix theory by a distribution of the inverse matrix-element variance. It keeps the basis invariance 
of the standard theory but violates the independence of the matrix elements. We consider the sub- 
extensive regime of q more than unity in which the transition from the Wigner to the Poisson 
statistics is expected to start. We calculate the level density for different values of the entropic 
index. Our results are consistent with an analogous calculation by Tsallis and collaborators. We 
calculate the spacing distribution for mixed systems with and without time-reversal symmetry. 
Comparing the result of calculation to a numerical experiment shows that the proposed non- 
extensive model provides a satisfactory description for the initial stage of the transition from chaos 
towards the Poisson statistics. 
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I. INTRODUCTION 



The past decade has witnessed a considerable interest devoted to non-conventional statis- 
tical mechanics. Much work in this direction followed the line initiated by Tsallis' seminal 
paper Q. The standard statistical —3 is based on the Shannon entropy _ 
S = —YtiPilnpi (we use Boltzmann's constant ks = 1), where {p{\ denotes the probabili- 
ties of the microscopic configurations. This entropy is extensive. For a composite system 
A + B, constituted of two independent subsystems A and B such that the probability 
p(A + B) = p(A)p(B), the entropy of the total S(A + B) = S(A) + S(B). Tsallis proposed 
a non-extensive generalization: S q = (1 — Eipf) /(q — 1). The entropic index q characterizes 
the degree of extensivity of the system. The entropy of the composite system A + B, the 
Tsallis' measure verifies 

S q (A + B) = S q (A) + S q (B) + (1 - q)S q (A)S q (B), (1) 

from which the denunciation non-extensive comes. Therefore, S q (A + B) < S q (A) + S q (B) 
if q > 1. This case is called sub-extensive. If q < 1, the system is in the super-extensive 
regime. The standard statistical mechanics recovered for q — 1. Applications of the Tsallis 
formalism eovered a w ld e elass of p_a; for a r ev le w please see, e.g. g However, 
the relation between the parameter q and the underlying microscopic dynamics is not fully 
understood yet. The value of q has been obtained from studies of dynamics in cases of 
ow- dimensional dissipative maps DO , and in some toy models of self-organized criticality 



Explicit expressions for q in terms of physical quantities exist in few cases, e.g. in 
turbulence problems 0] and physics of the solar plasma [jj. Aringazian and Mazhitov [^] 
obtained a Tsallis distribution function for a smaller subsystem weakly interacting with the 
remaining "quasi-thermostat" composed of a larger number M of particles, with an entropic 
index q — 1 ^ 1/M. 

A number of recent publications considered the possibility of a non-extensive general- 
ization to the random matrix theory (RMT) [9J. This is the statistical theory of random 
matrices H whose entries fluctuate as independent Gaussian random numbers. The matrix- 
element distribution has been obtained by extremizing Shannon's entropy subject to the 
constraint of normalization and existence of the expectation value of TrfH^H) What 
has become known as the Bohigas-Giannoni-Schmidt conjecture is that the quantum spectra 
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of classically chaotic systems are correlated according to RMT, whereas the spectral corre- 



lations of classically regular systems are close to Poissonian statistics [Tl|. Several attempts 
have been done to extend the applicability of RMT to include quantum systems with mixed 
regular-chaotic classical dynamics; for a review please see [l2j. For example, the principle of 
maximum entropy was used in for this purpose by introducing additional constraints con- 
cerning the off-diagonal elements ^j]. Non-extensive generalizations of RMT, on the other 
hand, extremize Tsallis' non-extensive entropy, rather than Shannon's. The first attempt in 
this d lr ect 10 n is p r oba Wy due to Evans and Michael Q. Toscano et a,. Q constructed 
non-Gaussian ensemble by minimizing Tsallis' entrop y an d obtained expressions for the level 
densities and spacing distributions. Bertuola et al. 16( have shown that Tsallis' statistics 
interpolates between RMT and an ensemble of Levy matrices 17] that have wide range 
of applications. They illustrated the spectral fluctuations in the sub-extensive regime by 
considering the gap function E(s) that gives the probability of finding an eigenvalue-free 
segment of length s. Analytical expressions for the level-spacing distributions of mixed sys- 
tems belonging to the three symmetry universality classes are obtained in jisj] . A slightly 
different application of non-extensive statistical mechanics to RMT is due to Nobre et al. 

Q. 

In this work, we use the integral representation of the gamma function to express the 
characteristics of the proposed non-extensive RMT in terms of integrals involving the char- 
acteristics of the conventional theory. We show that non-extensive statistics provides a 
principled way to accommodate systems with mixed regular-chaotic dynamics. 



II. NON-EXTENSIVE GENERALIZATION OF RMT 

RMT replaces the Hamiltonian of the system by an ensemble of Hamiltonians whose 
matrix elements are independent random numbers. Dyson showed that there are three 
generic ensembles of random matrices, defined in terms of the symmetry properties of the 
Hamiltonian. Time-reversal-invariant quantum system are represented by a Gaussian or- 
thogonal ensemble (GOE) of random matrices when the system has rotational symmetry 
and by a Gaussian symplectic ensemble (GSE) otherwise. Chaotic systems without time re- 
versal invariance are represented by the Gaussian unitary ensemble (GUE). The dimension 
/3 of the underlying parameter space is used to label these three ensembles: for GOE, GUE 
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and GSE, (3 takes the values 1, 2 and 4, respectively. Balian [h]] derived the weight functions 
Pp{H) for the three Gaussian ensembles from the maximum entropy principle postulating 
the existence of a second moment of the Hamiltonian. He applied the conventional Shannon 
definition for the entropy to ensembles of random matrices as S = — J dHPp(H) In Pp{H) 
and maximized it under the constraints of normalization of Pp(H) and fixed mean value of 
Tt(WH) . He obtained P P {H) oc exp [—r]Tr (H^H)~\ which is a Gaussian distribution with 
inverse variance l/2r/. In this section, we apply the maximum entropy principle, with Tsallis' 
entropy, to random-matrix ensembles belonging to the three canonical symmetry universali- 
ties. The Tsallis entropy is defined for the joint matrix-element probability density Pp{q, H) 
by 

S q [P fi (q, H)] =(l-JdH [Pp(q, I (q - 1). (2) 

We shall refer to the corresponding ensembles as the Tsallis orthogonal ensemble (TsOE), the 
Tsallis Unitary ensemble (TsUE), and the Tsallis symplectic ensemble (TsSE). For q — > 1, 
S q tends to Shannon's entropy, which yields the canonical Gaussian orthogonal, unitary or 
symplectic ensembles (GOE, GUE, GSE) 

There are more than one formulation of non-extensive statistics which mainly differ in 
the definition of the averaging. Some of them are discussed in j^jj. We apply the most 
recent formulation [22]. The probability distribution Pp(q,H) is obtained by maximizing 
the entropy under two conditions, 



/ 



dHP p {q ) H) = l 1 (3) 
<*- (4) 



JdH[Pp{q,H)] qr Ti(tfH) _ 2 
JdH[P,(q,H)f ^ 



where ap is a constant. The optimization of S q with these constraints yields a power-law 
type for Pp(q, H) 

Pp{q,H) = Z- 1 [1 + (q - 1% {Tr (H*H) - aj}]-^ , (5) 
where rj q > is related to the Lagrange multiplier rj associated with the constraint in (4) by 

V q = v/ J dH[Pp(q,H)}\ (6) 

and 

Z q = J dH [I + (q - l)rj q {Tr (H*H) - a 2 p }] ^ . (7) 



It turns out that the distribution (5) can be written hiding the presence of a\ in a more 
convenient form 

Pp(q, H) = Z' 1 [l + (q- l) Vq Tr (WH)] ~& , (8) 

where 



Vq JdH[Pp(q,H)F + (l- q ) V ay (9) 



and 



Z q = JdH[l + (q- l) Vq Tr (H^H)} . (10) 

The non-extensive distribution (8) is reduced to the statistical weight of the Gaussian en- 
semble when 5 = 1. 

It is important to note that the non-extensive distribution Pp(q, r) q , H) is isotropic in the 
Hilbert space because the dependence on the matrix elements of H enters through Tr (H^H) . 
In this way, Tsallis' statistics offers a random-matrix model for mixed systems, which is 
invariant under change of basis unlike most of the models in the literature. However, the 
distribution does not factorize into a product of distributions corresponding to the individual 
matrix elements if q ^ 1. Physically, this implies that the starting hypothesis of the standard 
RMT that the matrix elements are independent random variable does not hold in the non- 
extensive context described by Eq. (2). 

n 

The formalism developed in this section was applied in Ref. [18|] to ensembles of 2 x 2 
matrices. The calculation of the spacing distribution showed different behavior depending 
on whether q is above or below 1. It is found that the sub-extensive regime of q > 1 
corresponds to the evolution of a mixed system towards a state of order described by the 
Poisson statistics. On the other hand, the spectrum in the super-extensive regime develops 
towards the picked-fence type, such as the one obtained by Berry and Tabor [23] for the 
two-dimensional harmonic oscillator with non-commensurate frequencies. 



A. Sub-extensive regime 

In this paper, we shall consider only the sub-extensive regime, where q > 1. We note 

that Tt(H^H) = J2f=i ( H ii } ) 2 + 2 E?=iEi>j (^ 7) ) 2 ' where a11 the four matrices #W 
with 7 = 0,1,2,3 are real and where is symmetric while with 7 = 1,2,3 are 
antisymmetric. We introduce the new coordinates y = {yi, • • ■ , y^}, where d = N + (3N(N — 
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l)/2, and yf stand for the square of the diagonal elements or twice the square of the non- 
diagonal elements, respectively, and express the integrals in hyperspherical coordinates. The 
normalization condition (3) yields 



POO 

Z q (rj q ) = 2-^- 1 )/ 2 O d / y^dy [l + (q - l) Vq y 2 } ' 

Jo 
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r _ i 

= 2 -PN(N-l)/2 nd \Tl 2 J_ (11) 

2[(q-l) Vq f 2 T(i)V (jij) 

provided that q < 1 + 2/d , otherwise the integral diverges. Here = 2it d l 2 /Y(d/2) is the 
area of a unit <i-dimensional hypersphere and T(z) is Euler's gamma function. Condition 
(4), which now reads 

/°° __q_ I f°° q_ 

y d+1 dy[l + (q-l)r) q y 2 ] W J y^dy [l + (q - l) Vq y 2 ] ^. (12) 

The latter yield the following relationship between r] q and ah 

Vq = a 2 (2 + d - dq) ■ (13) 

For a Gaussian ensemble, <r§ = 2dv 2 , where v 2 is the variance of each of the non-diagonal 
matrix elements (or each of their components), so that 771 = l/(4t> 2 ). Condition (4) thus 
imposes the following upper limit on q 

<7<1 + ^, (14) 

beyond which the non-extensive formalism is not applicable for random matrix ensembles. 
This condition prevents the evolution of a chaotic system towards a state of order from 
reaching its terminal stage of the Poisson fluctuation statistics, as we may see in |l8| for 
the case of iV = 2, and later in this paper for the general case. The upper limit in (14) 
is essentially the extensive limit (q — > 1 + 0) since RMT is meant essentially for large 
matrices (d — > 00). For example, the sub-extensive regime for a GOE of 20 x 20 matrices 
is associated with values of q in the narrow range of 1 < q < 1.1. In spite of this, a minor 
non-extensivity produces a considerable effect on the spectral statistics of a large system as 
demonstrated below. This is attributed to the existence of the non-trivial "thermodynamic 
limit" N(q — 1) = constant, as pointed out by Botet et al. 
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B. Integral representation 



In the sub-extensive regime, the non-extensive RMT can arise from the extensive one by 
allowing the variances of the matrix elements to fluctuate using a transformation suggested 
by Wilk and Wlodarczyk and Beck 26 1. From Euler's representation of the Gamma 
function |27|, T(x) = J °° t^e^dt, one can easily derive the following expression 

[1 + (? - l)<Tr (H^H)r^ = — r t^e-A^-^H^^dt, (15) 

r (ft) Jo 

which is possible if q > 1. We now change the integration variable into rj — (q — l)rj q t. 
The joint distribution of matrix elements of a Tsallis random- matrix ensemble can then be 
expressed as 

Pf,(q,r, q ,H) = f°° f n ( V )^p\P/s( V ,H)drj, (16) 

JO L q KVq) 

P,(r ] ,H) = Z^ T < H ' H ), (17) 



where 



with 

z l( ,) = / (Ure ^) = ^^ ( 18 , 

is the distribution function for a Gaussian random-matrix ensemble with fluctuating matrix- 
element inverse variance 1/2^, and f n (f]) is the probability density of the ^-distribution (the 
distribution of sum of squares of n normal variables with zero mean and unit variance), 

-, / n n/2 / 

1 I 71 \ „7»/2-l I n V 



™=Tm\&) n ^v^j' (19) 

with order n = 2/(q — 1) and mean value r\ q = nd j\la\{n — d.)} = n /[Av 2 {n — d.)] . There- 
fore, the generalized distribution function Pp(q, rj q , H) of non-extensive statistics is expressed 
in terms of the distribution function Pp{f], H) of the corresponding Gaussian random-matrix 
ensemble by averaging over T], provided that r\ has a x 2 distribution. 

As mentioned above, that the non-extensive Hamiltonian matrix-element distribution in 
Eqs. (5) and (16) is invariant under change of basis. The mean value of each matrix element 
Hi]f = 0. On the other hand, the mean value of the square of a matrix element 



1 + Sa n ,„ . . 9 n — d 
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which is equal to the corresponding quantity in the standard RMT when n — > oo, as ex- 
pected. The distribution does not factorize into a product of distributions of individual 
matrix element, or matrix-element components, as in the standard RMT. The relative dis- 
persion of the squares of the matrix elements 




vanishes only in the extensive limit of n — > oo. We note that, for a given n and fixed v, the 
degree of correlation of matrix element measured by the covariance of their squares increases 
with increasing the dimension of the ensemble. This agrees with the result of non-extensive 
thermostatistics for the partition function of a system with N subsystems, which strongly 
suggest that the factorization approximation fails when N is large |29|. 



C. Eigenvalue distribution 



We now calculate the joint probability density for the eigenvalues of the Hamiltonian H. 
With H = U~ 1 XU, where U is the global unitary group, we introduce the elements of the 
diagonal matrix of eigenvalues X = diag(xi, • • • , xn) of the eigenvalues and the independent 
elements of U as new variables. Then the volume element (4) has the form 



dH= \ A N (X)fdXdn(U), 



(22) 



where (X) = [] n>m (i n — x m ) is the Vandermonde determinant and dfi(U) the invari- 
ant Haar measure of the unitary group 0, 12 1. The probability density Pp(H) is taken 
to be invariant under arbitrary rotations in the matrix space, Pp(j],H) = Pp{r],U~ 1 HU'). 
Integrating over U yields the joint probability density of eigenvalues in the form 

Y ^ — ^ poo / m \ d/2 



x N )dr], (23) 



where Pj^\rj,xi, ■■ ■ ,x^) is the eigenvalue distribution of the corresponding Gaussian en- 
semble, which is given by 



Pp (v, x u 



.xn 



Cp\A N (X)fexp 



A' 



i=l 



(24) 



S 



where Cp is a normalization constant. Similar relations can be obtained for the statistics 
that can be obtained from Pp(r) q , ; • ,xn) by integration. 



The fc-point correlation function |9j, [12J measures the probability density of finding a level 
near each of the positions x±, ■ ■ ■ ,x k , the remaining levels not being observed. It is obtained 
by integrating the eigenvalue joint probability density (16) over N — k arguments 

/oo /*oo 
dx k+1 --- dxffPjpfax!, ••• ,x N ). (25) 
-oo J — oo 

Therefore, the non-extensive generalization of the k-point function of a Gaussian ensemble 



R B k(Vi x u ■ • ■ , %k) is given by 



— 



where 

m = n — d 

Noting that ^- = 2v 2 m, we can easily see that the main parameter is m, which is subject 
to the restriction 

< m < oo. (27) 

The lower limit follow from the normalization condition as well as the constraint of finite 
average matrix norm (4). The upper limit corresponds the the standard RMT. 



III. LEVEL DENSITY 



The main goal of RMT is to describe the fluctuations of the energy spectra. Before 
the study of the fluctuations can be attempted, one must make a separation between the 
local level fluctuation from the overall energy dependence of the level separation. The level 
density of the standard random matrix ensembles is not directly related to the physical level 
density of the investigated systems. Nevertheless, it is essential to the proper unfolding of 
the spectral fluctuation measures. For the ^-dimensional GOE, the level density normalized 
to 1 is given by Wigner's semi-circle law 

pi (oo, e) = - ^jN^l - r]e 2 /N. (28) 

Here e — x/v is the energy expressed in units of standard deviation of the majority of 
matrix elements. In the following, we derive a corresponding formula for the non-extensive 
generalization of GOE, to which we shall refer as the Tsallis orthogonal ensemble (TsOE). 
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The level density is obtained by integrating the joint probability density of eigenvalues 
over all variables except one, so that the level density of TsOE p q (m,x) = 
Therefore, using Eq. (28) for izfj into Eq. (26), we obtain 

f ^ (2miV) m / 2 T (sg) fm + 1 m 2mN\ 
p g (m,e) = —= — 7 — r — , \ £ \Fi\ , h 2, — , 29 

HqK ! v^Fr (f ) r (f + 2) 1 1 V 2 2 £ 2 y v ; 



where \Fi (a,b,z) is Kummer's confluent hypergeometric function |2j|. At e — 0, Eq. (27) 
yields 

pJm,Q) = -\ - \V - (30) 

F9V ' ; vr V mA r (£) v ; 

Using the asymptotic properties of T(z) at large z |2j|, we can show that the height of 

p q (m,e) at the origin is lower than the GOE level density, pi(oo,0) = 1/ yKy/jfj, since the 

ratio of the gamma functions in the right-hand side of Eq. (30) is approximately equal to 

y/m/2[l - l/(4m)] for m 3> 1. At small m, where the relation F(l + z) = zY{z) tells us 

that T (m/2) pa 2/m, the dependence of p q (m,0) on m is mainly given by the factor yfm. 

On the other hand, the asymptotic behavior of p q (m,e) at large \e\ is given by 

p q {m,e) ~ |e|~ m_ . (31) 

Decreasing m from very large values lowers the magnitude of p q (m,e) at the origin below 
the semi-circle form of the Gaussian ensemble and raises its values at the periphery. This 
effect is clearly demonstrated in the left panel Fig. 1. The behavior shown in this figure is 
similar to that of the results of calculations by Toscano et al. Q|. 

In order to perform the statistical analysis of level fluctuations of the energy levels, 
one must take into account that the level density and hence the level spacing are strongly 
dependent on the intrinsic energy. For this purpose, the investigated spectra are transformed 
into the so-called 'unfolded' spectra j3] for which the local mean spacing is 1. On the other 
hand, calculations using RMT are performed for levels near the origin, where the level 
density is nearly equal to a constant proportional to y/rj. The energy scale is so far defined 
by the standard deviation of matrix elements, as, e.g. in Eq. (20). It is more suitable to 
express the quantities having the dimension of energy in terms of the mean level spacing 
rather than standard deviation of matrix elements. For this purpose we replace the ratio 
nrj/2rj q in Eq. (19) by 77/770 and define rjo by the requirement that the mean level spacing 
is 1. We then obtain for the non-extensive generalization any statistic Rp of a Gaussian 
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ensemble 

m 

9 r(f) J \r,J " u 'v 

where 77 is now understood as the square of the level density. 

IV. NEAREST-NEIGHBOR-SPACING DISTRIBUTION 

The nearest-neighbor-spacing distribution (NNSD) is frequently used for the analysis of 
experimental spectra. Unfortunately, RMT does not provide a closed form expression for 
NNSD. A very good approximation for this distribution is given by the so-called Wigner 
surmise which is the exact spacing distribution for Gaussian ensembles of 2 x 2 matrices. 
In this section we shall assume this approximation. We substitute the Wigner surmise for 
GOE and GUE into Eq. (32) and obtain expressions for NNSD of the corresponding Tallis 
ensembles. 



(32) 



A. Systems invariant under time reversal 

Chaotic systems, whose Hamiltonians are invariant under time reversal, are modeled in 
RMT by GOE. For these ensembles, the Wigner surmise is 

P GOE (r],s) =r]se~^ s \ (33) 

where rj is obtained by requiring that P GOE has a mean spacing equal 1. Substituting 
(26) into (25), we obtain the following expression for the non-extensive (Tsallis) orthogonal 
ensemble 

P T -°V,s) = k r h 'n - < 34 > 

(1 + |i7is 2 ) 



where 

"r(i(m-l)) 



7T 

Vl = o 



r(|m) 

is obtained by requiring that J °° sP Ts (m, s)ds = 1. 



(35) 



In the limit of m — > oo, r] m ir/m and the non-extensive NNSD approaches the extensive 



one, which is peaked at s — y2/ir « 0.80. At the other limit where m = 1, beyond which 
the mean the spacing distribution diverges. If one also requires that the second moment is 
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finite, one increases the lower limit into m = 2, for which one obtains 770 = tt 2 /2 and the 
NNSD becomes 

P™*(2,s)= f* 2 (36) 
(1 + n 2 s z ) 

which vanishes at the origin, has a peak at s — 2y / 2/ (tta/3) ~ 0.52, and decays asymptot- 
ically as s -3 . We therefore see that increasing non-extensivity modifies the NNSD from a 
Wigner form towards a Poisson distribution e~ s but never reaches it. Equation (34) agrees 
with the corresponding result obtained for 2x2 random matrix ensemble directly by in- 
tegrating the joint eigenvalue distribution This behavior is explicitly demonstrated 
in that paper. We thus expect the non-extensive RMT to describe the transition out of 
chaos, at least in its initial stage. We note that our result for the NNSD agrees with the 
corresponding distributions obtained in Ref. for the case of orthogonal universality. 

Figure 2 compares the NNSD in Eq. (34) with the corresponding results of a numeri- 
cal experiment 31]. This experiment imitates a one-parameter (denoted by S) transition 
between an ensemble of diagonal matrices with independently and uniformly distributed 
elements and a circular orthogonal or unitary ensemble. The figure shows that the TsOE 
distributions are generally in agreement with the numerical-experimental distribution al- 
though the quality of agreement gradually deteriorates, as expected, as the departure from 
chaos increases. 



B. Systems without time reversal symmetry 

RMT models systems whose Hamiltonians violates time reversibility, by GUE. The cor- 
responding Wigner surmise is 

pGUE/ n _ /2 3/2 J _W» 



(V) = ^/Ve-^ , (37) 

where 77 is obtained by requiring that P GUE has a mean spacing equal 1. Substituting (37) 
into (32), we obtain the following expression for the non-extensive (Tsallis) unitary ensemble 



^ £ Q (3 + m)) 4 /2 s* 

* r O) (1 + \V2S 2 ) 

where 



= V , 12 : :i , (38) 



r (m-1) 
7T T im 
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is obtained by requiring that J °° sP TsUE (m, s)ds = 1. The evolution of the distribution 
P TsUE (s) as m decreases from oo, where it is given by the Wigner surmise, a the limiting 
value of m = 2 is demonstrated in Ref. for the two dimensional case. Fig. 3 compares 
the NNSD for the TsUE with the corresponding distributions in the numerical experiment in 



Ref. 



31| . We again see the proposed non-extensive RMT provides a satisfactory description 



for the stochastic transition in terms of a single parameter, particularly in its early stage. 



V. CONCLUSION 



In the present work we have obtained a non-extensive generalization of the matrix-element 
theory by extremizing Tsallis's entropy, indexed by q, subjected to two constraints: normal- 
ization and finite average norm of the matrices. We consider the sub-extensive regime of 
q > 1, where the transition from chaos to order described by the Poisson statistics is ex- 
pected. The constraint of finite matrix-norm forces an upper limit of the entropic index:, 
limiting the attainable range tol<q<l + 2/d, where d is the dimension of the matrix- 
element space. This is essentially the extensive limit since RMT normally involves large 
matrices. Nonetheless, the obtained fluctuation statistics depend mainly on the parameter 
m = —d+2/(q — l), so that a minor deviation from extensivity leads to an observable effect. 
Because of this limitation, we expect the non-extensive formalism to provide a description 
for the early stages of transition from chaos towards regularity. We obtain distribution 
functions for the three symmetry universality classes, for which the probability of larger 
matrix elements decay algebraically instead of exponentially. By means of an integral trans- 
form, which is based on an integral representation of the gamma function, we express the 
characteristics of the non-extensive theory to those of the standard RMT. We calculate the 
level density and the NNSD for systems with mixed regular chaotic dynamics. We have 
also derived a generalization for the Wigner surmise that can be compared to numerical 
experiments with mixed systems. 
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Fig. 1. The level density (normalized to 1) for TsOE ploted against energy eigenvalues 
expressed in units of v for different values of the parameter m. 
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Fig. 2. NNSD for TsOE, calculated using Eq. 
numerical experiment by K. Zyczkowski and M Kus 
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Fig. 3. NNSD for TsUE, calculated using Eq. 
numerical experiment by K. Zyczkowski and M Kus 



'38) compared with the resuts of the 
3. 



18 




6 



